On the finite size behavior of quantum collective spin systems 



Giuseppe Liberti/'B Franco Piperno, 1,2 and Francesco Plastina 1,2 

1 Dipartimento di Fisica, Universitd della Calabria, 87036 Arcavacata di Rende (CS) Italy 
2 INFN - Gruppo collegato di Cosenza, 87036 Arcavacata di Rende (CS) Raly 

(Dated: July 5, 2009) 

We discuss the finite size behavior of the adiabatic Dicke model, describing the collective coupling 
of a set of iV-two level atoms (qubits) to a faster (electromagnetic) oscillator mode. The energy 
eigen-states of this system are shown to be directly related to those of another widely studied 
collective spin model, the uniaxial one. By employing an approximate continuum approach, we 
obtain a complete characterization of the properties of the latter, which we then use to evaluate 
the scaling properties of various observables for the original Dicke model near its quantum phase 
transition. 
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I. INTRODUCTION 

The interaction of N two-level systems (qubits) with a 
common single-mode quantum bosonic field is a paradig- 
matic example of collective quantum behavior. Dating 
back to the model put forward by Dicke, [l[, this has 
become one of the most investigated problems in quan- 
tum optics and condensed matter physics, with proposed 
physical implementations ranging from superconducting 
nano-devices, [H, to ultracold atoms and Bose- Einstein 
condensates in cavity Q. The Dicke model exhibits a 
second-order phase transition and, due to its broad 
application range [yL it has been studied extensively in 
the past few years, [1,0, HI, 0- It displays a rich dynam- 
ics, with many non-classical features (lCl UA L 12 . Hal ; in 
particular, the ground state entanglement [l4 15| and 
the Berry phase [TtJ of the Dicke Model have been 
diffusely analyzed and many aspects of its finite size be- 
havior have been obtained pH fl9L |2C| . 

The continued interest in the Dicke model also stems 
from the fact that it pertains to the same universality 
class as other intensely studied many-body systems that 
possess infinite-range interactions, and for which theo- 
retical models typically allow for exact solutions in the 
thermodynamic limit. 

The most general collective model of (effective) spin 
1/2 systems, the biaxial model in arbitrary field, can be 
described by the Hamiltonian (see Ref. [21( for details) 



k—x,y,z 



(i) 



where the Sk — J2i=i are collective Pauli op- 
erators that obey angular momentum-like commutation 
relations [Si,Sj] = 2ieijkSk- The energy eigenstates can 
be written in the angular momentum basis (we employ 
the standard one, apart from a factor 2 in the definitions) 
{|s, s z ); s z = — s, — s + 2, . . . , s — 2, s} constructed as the 



set of common eigenstates of both S 2 — S 2 + S 2 + S 2 

and S z . For a ferromagnetic interaction g x>y < 0, the 
ground state of the Hamiltonian belongs to the symmet- 
ric sub-space with S 2 — N(N + 2) and special and dif- 
fusely studied cases are the biaxial model in a transverse 
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0) Hjr Y (the well-known LGM model 

[H, HH) and the uniaxial model (S y = g y = 0) "• 

Under the thermodynamic limit, the phase diagram of 
these collective spin models has been simply established 
by a mean field approach (24J. For N large but finite, 
purely quantum effects become important and numeri- 
cal analysis have been implemented using the continuous 
unitary transformation method [25| and a semiclassical 
approach [2(i| . A qualitative understanding of the LGM- 
model is obtained in Ref. [13] introducing a double well 
structure above the phase transition in a semiclassical 
treatment of the system. 

In the present work, we establish an exact relationship 
between the Dicke model in the adiabatic regime (i.e. 
for the case of slow qubits coupled to a faster oscillator 
mode) and the uniaxial model, which is valid not only 
in the thermodynamic limit, but also for any finite num- 
ber N of spins. We then present an alternative analytic 
method which relies on a continuum approach to solve 
the collective uniaxial model for large N. We show that 
this method is useful to determine the finite size behavior 
and the entire 1/N expansion (i.e., critical exponents and 
pre-factors) at the critical point for both of the Dicke and 
the collective uniaxial spin models. These results corrob- 
orate several studies in which the exponents have already 
been derived. 

The objective of the present study is thus threefold, 
and the paper is organized accordingly: first, we consider 
the Dicke model in the regime in which the frequency of 
the quantum field is much larger than the energy spacing 
of the qubits; in this case, the field degree of freedom can 
be adiabatically separated from the qubit ones and an 
effective A^-qubit interaction can be obtained by means 
of the Born-Oppcnhcimcr approximation. This is done 
in SecJITl where the relationship with the uniaxial model 
is established for any energy eigen-state. Afterwards, we 
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focus on the quantum phase transition of this collective 
model (Sec. IIIII) for which we derive the 1/N expan- 
sion for some relevant physical observables and we also 
compute exactly various entanglement measures for the 
qubits. Finally, using these results together with those 
obtained in Sec. HTJ we obtain analogous 1/N expansions 
for the Dicke model (Sec. lIVp . A summary and some 
concluding remarks are finally given in Sec. [VJ 

II. ADIABATIC DICKE MODEL 

We consider a system of N qubits interacting with a 
single harmonic oscillator mode, described by the Hamil- 
tonian (h = c = 1) 

H=- 5 -S x + ^S z + utf a + a)S z (2) 
2 2 y A 

where a is the annihilation operator for the field mode of 
frequency co, S is the transition frequency of the qubit, e is 
the level asymmetry and A is the strength of the coupling 
between the oscillator and the two-level systems. 

We assume a slow qubit and work in the regime 
uj S> S by employing the Born-Oppenheimer approxima- 
tion, [3 [28| . The standard procedure is to separate the 
Hamiltonian of Eq.@ in t wo p arts, containing slow and 
fast variables, respectively [2!| 

H = H S + H } (3) 

where 

A S 
H f = Ljtfa+^S z + -=(tf + a)S z , H s = --S x (4) 
2 y N 2 

The eigenstates of the composite system can be written as 
a coherent superposition of the eigenkets of Hf, having 
a parametric dependence on (i.e. conditioned by) the 
values of the slow variables: 

IVV>= E ^ r VW). (r = 0,...,*) (5) 

{n,s*} 

where the displaced number states of the oscillator are 
given by 

\n[s z ]) =e-^ a *- a)s >\n}®\ S , Sz ), (6) 
They are the eigenstates of the fast Hamiltonian 

H f \n[s z }) =V n (s z )\n[s z ]), (7) 
with eigenvalues 

e A 2 

V n (s z ) = um+-s z - — s 2 . (8) 

For different n, V n (s z ) contribute an effective adiabatic 
potential felt by the slow subsystem so that the wave 
function </>i™' r ' ) of the A qubit system is determined by 

H eff ^ = E [ntr) (s z )^f\ (9) 



where the effective Hamiltonian is reduced to the form 

x 

Heff = ^S x + V n (S z ) = tvcJa + H%> 11 , (10) 

with H^' " being the Hamiltonian of the uniaxial model 
introduced in the previous section, with coupling con- 
stant g z = -j^: 

H^=- 6 -S x + ^S z + g z S 2 z . (11) 

The ground state of the coupled qubit-oscillator system 
is given by 

N 

\M= E ^ m e^ {af - a) \Q)<E>\N,m) (12) 

m=-N 

where tp m = 0m' O) . 

The uniaxial (as well as the LGM model) and the Dicke 
model are known to be equivalent in the thermodynamic 
limit. From the discussion of this section, we see that 
there is a strict relationship between the ground states of 
the two model-systems as both can be expressed in the 
angular momentum basis with the same amplitudes ip m . 
However, this last equation shows that at finite size there 
can be differences between their behaviors since, in the 
case of the Dicke model these coefficients gets effectively 
modified due to the presence of the displacement opera- 
tor, whose argument depends explicitly on the number of 
qubits N. This implies that we will find small differences 
in the 1/N expansions for the two models. 

In order to continue the discussion on the Dicke model, 
we need to evaluate the amplitudes ip m . Therefore, we 
now turn our attention to the uniaxial model with N 
qubits. Once the coefficients tp m are obtained, we will 
use them in Sec. IIVI to complete the description of the 
finite size behavior the Dicke model. 



III. UNIAXIAL MODEL 

A. Continuum approach 

The Hamiltonian of a uniaxial model for a spin system 
with a collective coupling can be written as 

H^=~S x + e -S z -^Sl (13) 

with 5 > and where we have re-scaled the ferromagnetic 
coupling constant by the number of spins, g z = —g/N. 

This is equivalent to the more diffusely found H^' " 
Hamiltonian, that can be obtained after the rotation 

To connect this model to the discussion of the previous 
section, one simply has to take g — A 2 /lu. 

The ground state of ' " lies in the maximum spin 
sector s = N. In this subspace, spanned by the states 
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{\N, to); to = -N, -N + 2, . . . ,N - 2,N}, the ground 
state can be written as 



iV 



I0O> = ^2 <Pm\ N , 



m) 



(14) 



where tp m are real coefficients. 

We limit our present discussion to the case of the sym- 
metric phase e = 0, that is the one relevant for the de- 
scription of the Dicke phase transition occurring at e = 
with A 2 = uS/4, corresponding to g — 6/4. 

In the angular momentum basis, the Hamiltonian takes 
a (N+1)®(N+1) tri-diagonal symmetric (Jacobi) matrix 
form with double symmetries along both the main and 
the second diagonal: 



T N = 



( X-N 





5-n 

X-N+2 
3-N+2 





&-N+2 
X-N+4 



\ 



where 



and 



5n-4 

$N-4 <5jV-2 
(5at_2 Ajv / 

(15) 



\ 2 

\ A 2 



-N < to < N 



■y/N{N + 2) - m(m + 2) . 



(16) 



(17) 



These coefficients satisfy a confluence property X m /N — > 
X(z ) and S 2 JN 2 — > 6 2 (z) where z = m/N as to, N — > oo 
[30| • By application of theorems on the zeros of orthogo- 
nal polynomials [3l| one finds that the ground-state en- 
ergy density in the to, ./V — > oo limit is given in general 
by 



£ (oo) = inf{A(z) - 2S(z)}. 



(18) 



Introducing the dimensionless parameter a = 4g/5 the 
minimum is found at 



0. 



(a<l) 



±^l~l/a 2 , (o>l), 



(19) 



and the corresponding thermodynamic limit of the 
ground state energy per spin is 



lim 



So(N) 

N 



-I, («<1), 
-£(« + £), (a>l). 



(20) 



For finite AT, the solution of the eigenvalues problem 
for the ground state reduces to the recurrence relation 



S m -2 l fim-2 + X m ip m + 5 m lfi 



m+2 



£0¥>r. 



(21) 



that can be rewritten as a second order linear difference 
equation 



2(5 m + 8 m -i) A 2 (fm + 2(d m 
+ (S m + <5 m -2 + X m )ip m = e Q ip ri 



Jm-2 



) Ai if,-, 



(22) 



where A 2 <^ m = {<p m +2 + fm-2 - 2tp m )/4 and Anp m = 
{ l Pm+2~ c Pm—2)/4 are finite differences of second and first 
order, respectively. 

A simple analytic behavior of the coefficients cp m for 
N ^> 1 can be derived by considering m/N as a contin- 
uous variable, and by expanding the recursion relation 
(f2"2")l in series around the minima of Eq. (|18p . For a < 1 , 
expanding in series (|2ip around to = and neglecting 
corrections of order 1 /N 2 , one obtains 



<Pr, 



e (N) 



1 



NS 2 
whose solution is 



1 



1 

N 



1 



4N 2 



'■Pr, 



2k \ 
ttNJ 



1/4 



-km 2 /4N 



(23) 



(24) 



with k — y/1 — a. The ground state energy per spin is 
given by 



£o(A0 



1 



1-VT= 



N 2 V N 

For a > 1, by expanding in series Eq. 
i ~ ±too = ±Ny/l — l/a 2 , one gets 



(25) 
IT) around 



a 
a 

27V 



1 

N8~ + 4 
a(a 2 — 1 
AN 2 



(to ± TOq)^ 



<p m ~0 (26) 



whose approximate solution is the symmetric superposi- 
tion 



'■Pr, 



T2 { ^ 



with 



2k 
^N 



1/4 



■<Pr, 



-fe(m=Fm ) /4N 



where k = a\J a 2 — 1. In this regime one has 
eo(N) < 



N 



a - 



N 



(27) 



(28) 



(29) 



In this language, the transition is readily understood: 
above the coupling value corresponding to a — 1 a dras- 
tic change in the form of the ground state wave function 
takes place, with a breaking of the "inversion" symme- 
try around to = 0. For a finite size system, the transition 
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0.20 




0.10 - 



200 



FIG. 1: (Color online): Normalized tp m function for a = 
0.3 (small red circles) and a — 1.3 (black circles) for a spin 
system of size N = 200. Comparison is made with the analytic 
expressions of Eqs, (l24l27|l (continuous lines). 



becomes smother and smother and the wave function cp m 
gradually changes from a one peaked gaussian to the su- 
perposition with two peaks that emerge progressively as 
the value of zq moves away from the origin (i.e., as a 
increases) . 

For large enough N, we can check the continuum ap- 
proximation by comparing it to the behavior obtained by 
solving the tridiagonal matrix numerically. In Fig. ([1]) (p m 
for N — 200 is shown with a — 0.3 and a = 1.3 compared 
with the analytic expressions of Eqs. (|24fl27|) . 



Finite-size corrections 



Having obtained the ground state coefficients ip m (to- 
gether with the ground state energy), we may evaluate 
the average values of every physical observable; in partic- 
ular, we concentrate on the total spin components. One 
immediately gets 



N 



2 de (N) 
' N dS 



1 



+ w( 1 + ^rfe)' 



NVa 2 -l ' 



(a > 1). 



(30) 



and 



(Si 
N 2 



1 de (N) 
' N dg 



(a < i); 



(31) 



The expression for (S 2 ) are, instead, a bit more com- 
plicated 



,o2 



= -[N(N + 2)-(S 2 z ) 



± 2 



N-2 

E 

m=-N+2 



at l a m <p m -2ip m +2 ■ (32) 



However, they can be simplified by making use of the 
simple results ^> m ~2^m+2 = e ~ 2k ^ N( Pm for a < 1 and 
<pm-2<Pm+2 = e~ 2k / N <pf n for a > 1 that are easily derived 
from our analytic expressions for tp m . Thus, one obtains 



(S 2 X 



1 + 77(1 



N 2 



N ax/a 2 -! : 



a < 1; 
a > 1. 



(33) 



(34) 



In the region a ~ 1 we must take into account also the 
next to leading order in the expansion of the recursion 
relation (f2Tj) that gives a non negligible contribution near 
the phase transition point. We thus need to consider the 
quartic-oscillator-like equation 




<Pm+ 



£o(AQ 
N8 



1 

N 



a 



A 1 



AN 2 



16N 4 



0. 



(35) 

Using the approach presented in a previous work [1 91 ] , the 
equation ([35]) can be reduced to a single-parametric prob- 
lem with the help of Symanzik scaling procedure [35[. 
This is done, by re-casting the equation ([55]) into the 
equivalent form 



Vn + ( e o (C) - (n 2 - n A ) tp n ~ 



(36) 



where n = m(2N)~ 2 / 3 is a scaled variable. The only 
remaining scale parameter is then £ = (2iV) 2 / 3 (l — a), 
while the ground-state energy is rewritten as 



6o(N) 
N 



For £ ~ (that is, very close to the transition point), we 
can resort to perturbation theory and obtain the ground 
state energy as an expansion in powers of £, 



eo(0 



OC 

E 

n=0 



(38) 



It is easy to show that /3q — eo(0) ~ 1.06036 is the lowest 
eigenvalue of the pure quartic oscillator and j3\ = e' (0) ~ 
0.36203. 

Using these results to obtain an approximate expres- 
sion for the ground state energy and for the coefficients 
ip m , it is easy to derive the following leading nontrivial 
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of the N spins in two blocks of size L and (N — L) , re- 
spectively. Using the decomposition 



-5 - 



-10 



-15 - 



<S. 2 >/N 2 



<S V 2 >/N 2 



6 

log 2 (N) 



10 



12 



FIG. 2: Scaling of two-spins correlation functions {S 2 )/N 2 
and (Sy)/N 2 as a function of N (log — log plot) at the critical 
point a = 1. 



finite-size corrections for one- and two-spin correlation 
functions 



N 



2/3i 



(2A0 2 /3 



(39) 



N,m) = p]^W-L,m-l)® \L,l) (43) 



l = -L 



where 



Pit: 



L 

L+l 
2 



N -L 

N-L+m-l 



N 

N+m 



(44) 



we obtain the ground-state reduced density matrix of the 
block of size L out of the total N spins in the form 



= E E \ L >h)(L,l 2 \ 

l 1 = -Ll 2 = -L 
N 

El/2 1/2 ,._>. 

P; 1 mP; 2 m-il+/2 < ^"^"^ i i+ i 2 ■ ( 45 ) 



=-JV 



We then compute the linear entropy as a measure of the 
entanglement of the block of size L with the rest of the 
system. 



(S 2 , 



N 2 



4/3i 



N 2 ~ (2iV) 2 / 3 



4/3i 



(2N) 2 / 3 



(40) 



(41) 



tl = Vl 



1 - Tr (pi 



N 



(46) 



2^-1 



9 

where the pre-factor is chosen to be 77^ = ttt — - in order 
to bound tl to 1. 

In particular, for L = 1, the state of every single qubit 
is found to be 



S(3 



TV 2 ~ 3(27V) 4 / 3 



(42) 



The critical exponents in these expressions are in full 
agreement with those reported in Ref.[25j]. Our method 
allowed us to obtain also the pre-factors, that cannot 
be determined with typical scaling arguments and that 
are important to transfer these results to the case of the 
Dicke model. 

In Fig. ^ we make a comparison of the analytical re- 
sults for the leading nontrivial finite-size corrections with 
those obtained from a direct numerical solution at the 
critical point. One can see that the agreement is good 
even for small values of N. 



1 ( r , (S x ) 



(47) 



where / is the identity. We can then evaluate the one- 
tangle as 



n 



1 - Tr (pl N ) 



= 1 - 



N 2 



One has 




a 2 JVdVo 2 -! : 



a < 1; 
a > 1. 



and 



(48) 



(49) 



C. Ground state Entanglement 

Before going back to the Dicke model, we use the re- 
sults we have obtained in order to discuss the critical 
behavior of the ground state entanglement for the uniax- 
ial model. In this respect, it is useful to make a partition 



n 



4/?i 



(2iV) 2 /3 



a = 1 



(50) 



The reduced density matrix of two qubits (L = 2), can be 
written in the angular momentum basis {|2, m) }, with 
to = 2, 0, —2. In general, one should also consider the 
state |0, 0); but its population is zero in our case, so that 



6 



we can erase the corresponding line and row, and write 
Pi n m the form: 



1.0 r 



P2,N = I y/2x^ 




(51) 



where the matrix elements may be expressed in terms of 
the expectation values of the collective operators as [12] 



v± 



± 



N(N - 2) + (S 2 ) ^ (S z ) 

AN(N-l) 
N 2 (S 2 ) 



2N 



4N(N- 



1 



N(N - 1) 



x± 



± 



[S+, S z ]- 



2N AN(N-l) 



(52) 
(53) 
(54) 
(55) 



The entanglement between two qubits can be expressed 
in terms of the concurrence [33| . Since the ground state 
lies in the maximum spin sector and has real coefficients 
in the basis {\N, m)}, one has 



C = max{0, C y } 



where 



(N-l)Cy = l- 



N 



(56) 



(57) 



Thus, the concurrence needs to be re-scaled by the factor 
N — 1 (thus, C vanishes ~ 1/N in the thermodynamic 
limit), with C r = (N — 1)C. In the thermodynamic limit, 
only C r remains finite: 

— — n-_ n- < 1 : 

(58) 

For finite N and at the critical point, Eq. ([42|) gives 

4/3 




C r d 1 



3(2JV)V3 



(59) 



This shows that, at the critical point, the behavior of 
the concurrence is modified and C r scales with N with a 
critical exponent of 1/3. 



IV. FINITE SIZE SCALING OF THE DICKE 
MODEL 

We can now make use of the expressions of the ampli- 
tudes ip m obtained for uniaxial model in order to discuss 
the finite size behavior of the Dicke model. 

Whenever we are interested in a qubit observable, 
that is, whenever the result can be obtained by trac- 
ing out the oscillator, the only difference with the uniax- 
ial model is the appearance of exponentials of the form 




FIG. 3: (Color online): Scaled concurrence for the Dicke 
model as a function of a, for D — 0.1 and for system sizes 
N = 10, 20, 40, 100 and oo (bottom to top). 



exp{— A 2 (to - 



S\2 



/Nlo 2 }, due to the overlap of different 



coherent states. This kind of terms modifies the behavior 
of the spin observables for small N, but for very large N 
one can expect to obtain very similar behaviors for the 
Dicke and the uniaxial models. One finds 



and 



(S 2 S 2 } — > e 



(S x 



■(sl 



(60) 



(61) 



which coincides, respectively, with Eq. (|30| and Eq. (J32J) 
in the D = S/uj — > limit. Using these results, one can 
show, for example, that the ground state energy at the 
critical point reads 



5 , 

Y "-2 11 



D 



Wo 



2N 



(2iV) 4 /3 



(62) 



Once all of the average values of the spin observables 
are obtained, it is easy to get expressions for the vari- 
ous entanglement measures. In particular, the re-scaled 
concurrence in the thermodynamic limit, reads 



(63) 




where a = (1 + D 2 )/2D = (S 2 - 
size, at the critical point one gets 



u) 2 )/2uj5. For finite 



1 — D 



4/3 



3(2iV)i/3 



(64) 



Thus, the concurrence scales with N exactly as in the 
uniaxial model. Fig. shows the C r both for finite N 
and for N — > oo. 

The difference between the adiabatic Dicke model and 
the uniaxial one lies in the presence of the oscillator, 
which is far detuned from the spins but that can still be 
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FIG. 4: (Color online): The tangle rjv between the oscillator 
and the N qubits as a function of a, for D = 0.1 and for 
system sizes N = 6, 10, 20, 100 and oo (bottom to top). 



excited (because of the presence of the counter-rotating 
terms in the Hamiltonian) and becomes correlated with 
the qubits. In particular, the entanglement between the 
oscillator and the N qubits, can be evaluated by the lin- 
ear entropy which is of the form 



tn = r] N (l-Tr{p N }^ 



(65) 



where pm is the reduced density matrix for the N-qubits 
sub-system, obtained from the ground state density op- 
erator (|12p by tracing out the field variables 

Pn = TM |^o>(Vo|} (66) 

N 

J2 e-^ m ^ m ^^ mi ^ n2 \N 7 m 1 )(N 7 m 2 \ 

mi 7 m2 — — N 

Evaluating the trace of pn squared, one has 



TN 



T) 1 



N 

E 

mi ,7712= — N 



e-^ (mi - l2)2 <^ 2 | (67) 



In the thermodynamic limit the sum can be computed 
exactly to get 



1-1 



1 




(a<l) 
(a>l), 



(68) 



which shows a cusp at the critical point, where = 1. 
Fig. (HJ) shows tjv both for finite N and for N — > oo. 
When N is very large (N ^> 4/D 3 ), the entanglement 
scales as 



T N (a = 1) 



1/6 



(69) 



where K — \ J dnip^ ~ 0.46, and <p n is the normalized 
solution of Eq. for a = 1. The fact that the lead- 
ing term in the 1/N expansion of tn has exponent 1/6, 



implies that the convergence of the series is slower (with 
respect to those found for other physical quantities) and 
that for small values of N, subsequent terms should be 
taken into account. 



V. CONCLUDING REMARKS 

We have discussed the finite size critical behavior of 
the Dicke model for the case of a fast oscillator coupled 
to many slower qubits. We have derived a direct rela- 
tionship between this system and the uniaxial model, de- 
scribing the collective interaction among qubits residing 
on a fully connected graph. In particular, we have ob- 
tained a precise one-to-one correspondence between the 
energy eigen-states of the two models, showing that their 
critical behavior are closely related both in the thermo- 
dynamic limit and at finite size. We have then, adopted 
a continuum approximation in order to describe analyt- 
ically the ground state of the uniaxial model which we 
used to re-obtain all the known features of the model, 
such as its critical exponents. We have also obtained a 
full characterization of the 1 /N expansion (including non 
universal features such as the pre-factors) for many phys- 
ical observables, among which we dedicated a particular 
emphasis to the description of the entanglement content 
of the ground state and to its critical behavior. 

Using the solution obtained for the uniaxial model, we 
were then able to go back to the original Dicke model and 
to describe its critical behavior and its scaling properties, 
again obtaining not only the scaling exponents for various 
physical quantities, but their entire 1/N expansions (of 
which the first terms are shown and discussed explicitly). 

The two models we have discussed obviously differ be- 
cause of the presence of the bosonic mode in the Dicke 
case. From a physical point of view this implies that an 
entanglement is built up not only among the qubits (as 
in the uniaxial model) but also between qubits and oscil- 
lator. Formally, this manifests itself in the fact that the 
oscillator state is a displaced vacuum (i.e. coherent) state 
conditioned on the qubit magnetization in the direction 
of the coupling (i.e., on the value of S z , in our notation). 
The presence of these quantum correlation with the oscil- 
lator also modifies the entanglement among qubits (for- 
mally, because of the presence of some exponential pre- 
factors that essentially suppress entanglement), and this 
can be interpreted in terms of the monogamy of entan- 
glement. 

Apart from this, the two models have many features in 
common; in particular, their critical behavior are closely 
related and their quantum phase transition are essen- 
tially the same, occurring at the same point in parameter 
space (once the proper relationship between the physical 
parameters is taken into account). 
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